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SUMMARY 

Non-linear programming algorithms play an important role in structural design optimization. Fortunately, 
several algorithms with computer codes are available. At NASA Lewis Research Centre, a project was 
initiated to assess the performance of eight different optimizers through the development of a computer code 
CometBoards. This paper summarizes the conclusions of that research. CometBoards was employed to 
solve sets of small, medium and large structural problems, using the eight different optimizers on a Cray- 
YMP8E/8128 computer. The reliability and efficiency of the optimizers were determined from the perfor- 
mance of these problems. For small problems, the performance of most of the optimizers could be 
considered adequate. For large problems, however, three optimizers (two sequential quadratic programming 
routines, DNCONG of IMSL and SQP of IDESIGN, along with Sequential Unconstrained Minimizations 
Technique SUMT) outperformed others. At optimum, most optimizers captured an identical number of 
active displacement and frequency constraints but the number of active stress constraints differed among the 
optimizers. This discrepancy can be attributed to singularity conditions in the optimization and the 
alleviation of this discrepancy can improve the efficiency of optimizers. 
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INTRODUCTION 

Non-linear programming algorithms play an important role in structural design optimization. 
Fortunately, several algorithms with computer codes have been developed during the past few 
decades. To assess the performance of different optimizers, a project was initiated at NASA Lewis 
Research Centre and a computer code called CometBoards, which is an acronym for Compara- 
tive Evaluation Test Bed of Optimization and Analysis Routines for the Design of Structures 1 has 
been developed. CometBoards incorporates eight popular optimization codes: the Sequential 
Unconstrained Minimizations Technique (SUMT); 2 the Sequential Linear Programming method 
(SLP); 3 the Feasible Directions method (FD); 3 the Sequential Quadratic Programming technique 
(SQP of IDESIGN); 4 the DNCONG of the IMSL routine; 5 NPSOL, which is available in the 
NAG library; 6 the Reduced Gradient method (RG); 7 and the Optimality Criteria methods (OC). 8 
CometBoards was employed to solve a set of 41 structural problems by using its optimizers on 
a Cray-YMP8E/8128 computer. The reliability and efficiency of the eight optimizers were 
ascertained on the basis of the performance of these problems. The problems were solved for 
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multiple load conditions, and behaviour constraints were imposed on stresses, displacements and 
frequencies. The examples were selected so that at optimum, numerous stress, displacement and 
frequency constraints were active. Initial design, upper and lower bounds and convergence 
parameters were specified to ensure that the evaluation had no bias towards any particular 
optimizer or any particular problem. The eight optimizers might have been updated during the 
time CometBoards was developed, but any such improvements were not accounted for. 

Evaluations of optimizers that are available in the literature 9 " 18 deal broadly with individual 
code validation by their developers. The studies lack uniformity because problems and computa- 
tional platforms differ and the evaluations are over a decade old. For example, Arora, the 
developer of SQP of IDESIGN, compared his algorithm to the NAG/NPSOL optimizer. 9 11 
Most of Arora’s problems were trusses for stress and displacement constraints and were solved on 
a PRIME 750 computer. Schittkowski, who is the developer of the DNCONG optimization 
routine in the IMSL library, essentially validated his code 12 " 14 by solving many theoretical 
examples. Venkayya, one of the developers of ASTROS in which OC and FD optimizers are 
used, 16 attempted an evaluation of a few practical problems on a VAX 11/785 computer. 15 An 
intermediate complexity wing problem, used by Venkayya with stress and displacement con- 
straints, 15 is also included in our test bed with the addition of frequency constraint. Ragsdell’s 
evaluation, 17,18 includes mostly simple mechanical application problems. The current paper 
differs from those available in the literature in several respects: (1) a single tool, CometBoards, 
evaluates all eight optimizers on a common Cray-YMP computer; (2) solutions to a set of 
problems, which were grouped into categories of small, medium and large, are used; and (3) design 
parameters were selected to ensure that the evaluation had no bias towards problems or optimizers. 
In brief, the comprehensive evaluation presented in this paper does not duplicate previous work. 
This paper presents a brief theory of optimization methods, a description of CometBoards, 
a summary of the numerical examples and their solutions, discussion and conclusions. 

THEORY OF OPTIMIZATION METHODS 

Structural design can be formulated as: Find the n design variables £ within the prescribed upper 
and lower bounds (x^ x, xF, i = 1,2,... ,n) which make a scalar objective function /(Y) an 
extremum (here, minimum weight) subject to: a set of mi inequality constraints gj( 3f) ^ 0, 
(j = 1 , 2 , . . . ,OTj) and m c equality constraints g J+mi (X) = 0 (j =1,2,... ,m t ). 

Stress, displacement and frequency behaviour constraints were considered in this study. 
A cursory account of the different optimization methods available in CometBoards is provided 
herein. Readers may refer to specified references for details. 

SUMT, as implemented in the code NEWSUMT, is available in CometBoards. In NEWS- 
UMT, the penalty function has been modified to improve efficiency and a modified Newton’s 
approach is used to calculate the direction vector while a golden section technique is used to 
determine step length. SLP, as implemented in the Design Optimization Tools (DOT 2.0), 
is available in CometBoards. From the original non-linear problem, a linear programming 
subproblem is obtained by linearizing a set of critical constraints and the objective function 
around a design point. The linearization process and linear solution sequence is repeated until 
convergence is achieved. FD, as implemented in DOT 2.0, is available in CometBoards. In FD, 
a direction is obtained that is both usable and feasible. A minimum along the search direction is 
generated by polynomial approximation. Three implementations of the sequential quadratic 
programming technique, SQP of IDESIGN, DNCONG of IMSL and NPSOL in NAG, are 
available in CometBoards. In this technique, the original non-linear problem is solved 
through a sequence of quadratic subproblems. In SQP of IDESIGN, a Lagrangian function is 
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approximated. The step length is obtained by minimizing a composite descent function. 
DNCONG of IMSL uses quasi-Newton updates for the Hessian of the Lagrangian function while 
the constraints are linearized. 12, 14,19 The step length for an augmented Lagrangian is calculated 
using a bisection method. 20 NPSOL in NAG also uses an augmented Lagrangian. The search 
direction is generated through a quadratic subproblem while step length is calculated using an 
augmented Lagrangian, which is designed to avoid discontinuities as much as possible. RG, as 
implemented in the code OPT, 7 21 " 23 has been incorporated into CometBoards. This method 
partitions the design variable into decision and slave variables and a reduced gradient is used to 
generate a search direction. A line search is carried out by bounding the minimum and then 
calculating the minimum within some tolerance. 24 ' 26 OC, available in CometBoards, can be 
considered as a variant of the Lagrange multiplier approach applied to structural design problems. 
In OC, an iterative scheme is followed to update the multipliers and the design variables separately. 

In addition to the eight popular optimizers (SUMT, SLP, FD, SQP in IMSL, SQP of 
IDESIGN, NPSOL in NAG, RG and OC) considered in this paper, there are other routines 
which can be used for design applications. These other routines include, mFD (modified FD), 27 
GRG2 (Generalized Reduced Gradient method), 28 CONLIN, 29 MOM (Method Of Multi- 
pliers), 30 Genetic, 31 Convex, 32 Asymptotes, 33 dual methods, 34 etc. Mention has been made 
regarding the potential of genetic and asymptote algorithms. However, to date, the merits for 
both of these algorithms have been shown for only rather modest problems. For example, the 
potential of the asymptote algorithm 33 is shown through solutions for a determinate cantilevered 
beam with a single displacement constraint and for two simple trusses with only a few bars for 
stress limitations. In brief, the maturity and applicability of some of these algorithms for complex 
structural systems, under static and dynamic constraints, still remain to be proven, and at this 
time these algorithms have not been incorporated into the CometBoards test bed. 


DESCRIPTION OF COMETBOARDS 

The basic organization of CometBoards is depicted in Figure 1. The central executive with 
command level interface (Figure 1) links the three modules (optimizer, analyser and data input) 


CometBoards 



Figure 1. Comparative evaluation test bed of optimization and analysis routines for the design of structures (CometBoards) 
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of the code to formulate and solve an optimization problem. Note that there are eight choices 
for optimizers. The analyser options are the displacement method, 8, 35 the integrated force 
method, 8,36 the simplified force method, 8 etc. There are three input data files, one for analysis 
(anldat), one for design (dsgndat) and one for optimization (optdat). CometBoards has consider- 
able flexibility in solving a design problem by choosing any one of the eight optimizers and any 
one of the three analysers. A more detailed description of CometBoards can be found in Reference 1. 


EXAMPLE PROBLEMS 

The numerical testbed of CometBoards includes over 41 problems, most of which were taken 
from the literature. 1,8, 15, 37-42 Minimum weight was the objective and a linking strategy was 
followed to reduce the number of design variables. Stress, displacement and frequency behaviour 
constraints were considered. Multiple static load conditions and consistent elemental mass for 
dynamic analyses were also considered. The load conditions, mass distributions, and behaviour 
limitations were selected to ensure that several types of behaviour constraints were active at the 
optimum. The initial design of unity was considered for all problems unless otherwise specified. 
A consistent set of upper and lower bounds was specified for each problem. Typically, default 
optimization parameters and convergence criteria specified in the individual codes were used. 
These parameters, however, were changed when convergence difficulty was encountered. Results 
for all 41 examples are summarized in Table I. The normalized optimum weight and the 
normalized Cray-YMP8E/8 128 CPU time for a select set of 14 examples are given in Table II and 
depicted in Plates 1 through 4. The weight was normalized with respect to the optimum weight 
obtained for the best feasible design. A brief description of the 14 examples follows. 


Examples (P\a-P\d). Three-bar truss 

The popular three-bar truss 8, 37, 41 (with modulus E = 30000 ksi and density p = 0-1 lb/in 3 .) 
was subjected to a single load condition. It had three design variables and six constraints (three 
stress, two displacement and one frequency). Optimum weight and CPU time are depicted in 
Table II (Pla, Plb, Pic, Pld) and Plates 1 and 4. The optimum weight was 92*87 lb and three 
constraints (one stress, one displacement and one frequency) were active. Seven optimizers 
(SUMT, SLP, FD, SQP, IMSL, NPSOL and RG) performed satisfactorily. OC was inadequate, 
yielding a 38-6 per cent overdesign. The problem was solved again for three different initial 
designs (the SUMT optimum design, 150 per cent of SUMT optimum and 50 per cent of SUMT 
optimum). Results followed the earlier pattern where the initial design was unity. The CPU times 
on the Cray-YMP computer required for different optimizers are depicted in Plate 4. For unit 
initial design, SLP required the least CPU time of 0 07 s, while RG was most expensive at 3*18s. 


Example P2 . Tapered 10-bar truss 

A tapered 10-bar aluminum truss 8 was subjected to two load conditions. It had 10 design 
variables and 25 behaviour constraints (20 stress, four displacement and one frequency). The 
optimum weight was 3326-74 lb with 1 1 active constraints (eight stress, two displacement, and one 
frequency). Four optimizers (SUMT, SQP, IMSL and NPSOL) converged for this example. SLP 
and FD converged to an 82-4 per cent under-design condition. Optimizer RG failed and OC was 
marginal at a 5-6 per cent over-design. Cray-YMP CPU time varied between 1*28 s for NPSOL 
and 1 91 s for SUMT. 



Normalized weight (1 = optimum) 



Plato I Performance of different optimizers for small problems 
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Plate 2 Performance of different optimizers for medium problems 
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Plate 3 Performance of different optimizers for large problems 
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Plate 4 Cray-YMP CPU time for different optimization methods for 14 problems 
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Example P3. Tapered cantilever beam 

The cantilever truss of Example P2 was modelled next using eight triangular membrane 
elements. 43 The loads and constraints were kept identical to those in Example P2. There are eight 
thicknesses of the elements which were considered the eight design variables. The problem had 21 
constraints (16 von Mises stress, four displacement and one frequency). Optimum results ob- 
tained are given in Table II. The optimum weight was 1440-24 lb with six active stress constraints. 
Seven optimizers (SUMT, SLP, FD, SQP, IMSL, NPSOL and RG) performed well while OC 
produced a 2-8 per cent over-design (see Table II). Cray-YMP CPU time varied from 1-62 s for 
SLP to 11 -22 s for RG. 


Example P4. 25-bar truss 

A 25-bar aluminum truss 37, 38 had eight linked design variables and was subjected to two load 
conditions. It had a total of 86 behaviour constraints (50 stress and 36 displacement). Six 
optimizers (SUMT, SLP, FD, SQP, IMSL and NPSOL) converged to an optimum weight of 
544-73 lb with four active displacement constraints (Tables I and II). Optimizers RG and OC 
failed. Cray-YMP CPU time ranged from 1-64 s for SQP to 615s for NPSOL. 


Example P5. 165-ft-tall antenna tower 

A 165-ft-tall steel antenna tower with 252 members 37 had six linked design variables and was 
subjected to two load conditions. Its overhead dish antenna was modelled as a lumped mass for 
frequency calculations. It had a total of 529 behaviour constraints (504 stress, 24 displacement 
and one frequency). Five optimizers (SLP, FD, SQP, IMSL and NPSOL) converged to an 
optimum solution of 5299-84 lb with small deviations (Table II). At the optimum, six stress, 12 
displacement and one frequency constraints are active. Optimizers RG and OC failed while 
SUMT produced a 6 per cent under-design. The Cray-YMP CPU time varied between 222-71 s 
for SLP and 1893-80 s for NPSOL. 


Example P6. 60-bar trussed ring 

A 60-bar trussed aluminum ring 8 was subjected to three load conditions and had two lumped 
masses. It had a total of 184 constraints (180 stress, three displacement and one frequency) and 25 
linked design variables. The optimum weight was 414.51 lbs, and at optimum, 22 stress, one 
displacement and one frequency constraints were active. Six optimizers (SUMT, SLP, FD, SQP, 
IMSL and NPSOL) converged (Table II). Optimizer RG failed, whereas OC produced a 4-1 
per cent over-design with a IT per cent constraint violation. Cray-YMP CPU solution time 
ranged from 12-67 s for SLP to 144-11 s for NPSOL. 


Example P7. Geodesic dome 

A geodesic dome, 39, 40 with a diameter of 240 in. and a height of 30 in., was subjected to a single 
load condition. It was modelled using 156 bars and 96 triangular membrane elements. The bars 
were made of a material with modulus E = 30000 ksi, and density p = 0-1 lb/in 3 . Membranes 
were made of aluminum, with modulus E = 10 000 ksi, and density p = 01 lb/in 3 . The bar areas 
and membrane thicknesses were grouped to obtain eight and four linked design variables, 
respectively. The dome had a total of 254 constraints (156 stresses for bars, 96 von Mises stresses 
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for membranes, one displacement and one frequency). The optimum weight obtained was 
1022*67 lb with 170 active constraints (168 stress constraints, one displacement and one frequency 
(see Table I). Four optimizers (SUMT, SQP, IMSL and NPSOL) converged with small devi- 
ations. Optimizers RG and OC failed, while SLP and FD produced over-designs (12*7 and 12*0 
per cent, respectively). The Cray-YMP CPU time varied between 448*32 s for SUMT and 548*36 s 
for NPSOL. 


Examples (P8a~P8d). Intermediate complexity wing 

An intermediate complexity wing 8, 15 was modelled with a total of 158 elements consisting of 
39 bars, two triangular membranes, 62 quadrilateral membranes and 55 shear panels. The wing is 
made of aluminum with modulus E = 10 500 ksi and density p = 0*1 lb/in 3 . The elements were 
grouped to obtain 57 linked design variables. The wing, which was subjected to two load 
conditions, had a total of 321 behaviour constraints (316 stress, four displacement and one 
frequency). The optimum design for this problem was obtained from four different initial points: 
(1) initial design of unity; (2) initial design equal to the SUMT optimum design; (3) initial design 
equal to 150 per cent of the SUMT optimum design; and (4) initial design which is infeasible at 50 
per cent lower than the SUMT optimum design. Results obtained for all four cases are 
summarized in Table II (P8a-P8d). The optimum design was 387*76 lb and there were a total of 
119 active constraints (117 stress, one displacement and one frequency). For initial design equal to 
unity (see Table II, Problem P8a), optimizers SUMT, FD, SQP, IMSL and NPSOL reached the 
optimum within a 3*7 per cent error margin. Optimizer RG failed to solve the problem. 
Optimizers SLP and OC also failed to converge to the optimum (producing 14*3 and 20*1 per cent 
over-designs, respectively). Cray-YMP CPU time varied between 1075*21 s for SQP and 8292*35 s 
for NPSOL. 


DISCUSSION 

For the purpose of this discussion, the 41 examples of the CometBoards test bed are grouped as 
small, medium and large problems. The number of linked design variables ranged between 3 and 
19 for small problems (Group I). Group I contains a total of 19 problems, which are designated as 
Pla-P5, P7 and P9-P18. The normalized optimum weight for small problems obtained by each 
optimizer is depicted in Plate 1. For medium problems (Group II), the number of linked design 
variables ranged between 20 and 39. There are 12 medium problems, which are designated as P6, 
P19-P23 and P30-P35. The normalized optimum weight for the medium problems obtained by 
each optimizer is illustrated in Plate 2. Problems with more than 40 independent design variables 
are referred to as large problems (Group III). There are 10 large problems which are designated as 
P8a-P8d, and P24-P29. The normalized optimum weight for large problems obtained by each 
optimizer is depicted in Plate 3. 

The discussion is separated into the following five categories: (1) convergence to the optimum 
weight, (2) number of active constraints at optimum, (3) Cray-YMP8E/8128 CPU time required 
to solve the problem, (4) singularity in structural optimization and (5) default optimization 
parameters. 


Convergence to optimum weight 

The normalized optimum weights for all 41 problems, obtained by the eight optimizers, are 
depicted in Plates 1-3 for small, medium and large problems, respectively. In these Plates, unity 
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represents optimum weight and more than unity indicates over-design, while less than unity is 
infeasible design. For the purpose of comparison, a solution with constraint violation of less than 
1 per cent and weight which is within 1 per cent of the best feasible design is considered correct. 
A design is acceptable when the constraint violation is less than 1 per cent and the weight is within 
5 0 per cent of the best obtained by the eight optimizers. Convergence to optimum weight for each 
of the eight optimizers follows. 

(1) SUMT converged to optimum solution for 35 of 41 examples, which consisted of 17 small, 
nine medium, and nine large problems. SUMT failed for four problems. These are: one 
small problem (P5), two medium problems (P21 and P23) and one large problem (P8b). For 
the two medium problems, the SUMT solution was more than 1 per cent infeasible. For the 
large problem, SUMT gave an over-design of more than 5 per cent. 

(2) SLP of DOT 2.0 successfully solved 19 of 41 examples, which consisted of 11 small, seven 
medium, and one large problems. SLP failed for seven small problems (P2, P7, P9, P10 and 
P16-P18), three medium problems (P19, P21 and P23) and nine large problems (P8a, P8c, 
P8d and P24-P29). 

(3) FD of DOT 2.0, successfully solved 16 of 41 examples, which consisted of nine small, six 
medium and one large problems. FD produced infeasible designs for six small problems 
(P2, P9, P10 and P16-P18) and one large problem (P8d). It produced over-design condi- 
tions for eight problems (P7, P8c, P19, P24 and P26-P29). 

(4) SQP of IDESIGN successfully solved 32 of 41 examples, which consisted of 15 small, 10 
medium and seven large problems. This optimizer failed to give a feasible optimum design 
for three small problems (P15, P17 and P18), two medium problems (P19 and P22) and 
three large problems (P8c, P8d and P25). 

(5) IMSL optimizer DNCONG successfully solved 37 of 41 examples, which consisted of 17 
small, 12 medium and eight large problems. DNCONG of IMSL failed to optimize the 
intermediate complexity wing (Problems P8c and P8d). 

(6) NPSOL successfully solved 25 of 41 examples, which consisted of 13 small, eight medium 
and four large problems. This optimizer failed (with a 1 per cent infeasible constraint) for 
four small problems (P15-P18), four medium problems (P23, P31, P32 and P34) and four 
large problems (P8d, P25, P26 and P28). It produced more than 5 per cent over-design for 
large problem P8c. 

(7) RG successfully solved 13 of 41 examples, which consisted of seven small, four medium and 
two large problems. RG failed for 12 small problems. It also failed for seven medium 
problems and three large problems. The optimizer RG failed with well over 100 per cent 
error in the optimum weight for 15 problems. 

(8) OC successfully solved 16 of 41 examples, which consisted of six small, five medium and five 
large problems. OC failed for nine small, two medium and five large problems with an error 
in the optimum weight exceeding 5 per cent, as well as for three medium problems with an 
infeasible design greater than 1 per cent. 


Number of active constraints at optimum 

The number of active constraints at the optimum for all examples is given in Table I. Typically, 
different optimizers produced identical numbers of active frequency and active displacement 
constraints. However, the number of active stress constraints generated depended on the opti- 
mizer of choice. For example, with the geodesic dome problem (P7), the number of active stress 
constraints produced were 168 by SQP of IDESIGN, 162 by SUMT and NPSOL, and 156 by 
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IMSL. Consider next the set of five examples depicted in Table III that failed to converge. The 
minimum weights ranged between 3-2 and 12-7 per cent over-designs or under-designs. These 
examples produced correct numbers of displacement and frequency constraints, but failed to 
produce the correct numbers of active stress constraints. The deficiency in the number of active 
stress constramts ranged between 3 for Problem P2 and 42 for Problem P8a. For these problems 
the failure of the optimizers could be attributed to their inability to produce the correct number of 
active stress constraints. This aspect is also described in the section entitled, ‘Singularity in 
structural optimization’ of this paper. 


CPU time required for the solution 

The normalized CPU times on a Cray-YMP8E/8128 computer were recorded for a set of 14 
examples. The normalization was with respect to SQP of IDESIGN except for Problems P8c and 
P8d, which were normalized with respect to SUMT (Table II and Plate 4). CPU time required to 
solve the same problem differed among optimizers. Even for a small problem (Pla), normalized 
CPU time differed from 0 507 for SLP to 22 069 for RG. For a medium problem (P6), normalized 
time differed between 0-343 for SLP to 3-899 for NPSOL. For a large problem (P8a), normalized 
CPU time varied from 1-695 for IMSL to 1-116 for SUMT and 1 000 for SQP of IDESIGN. It is 
observed that variation in CPU time was rather mild for large problems. 


Singularity in structural optimization 

Singularity in structural optimization was identified for three situations: 8,41 (1) The number of 
active constraints exceeds the number of design variables. Out of the 41 test bed problems, the 14 
examples listed in Table IV are prone to this type of singularity. (2) Linear functional depend- 
encies exist among a small number of active stress constraints. This type of singularity is 
suspected to have occurred for some of the examples given in Table III. (3) Linear functional 
dependencies exist among a small number of active stress and displacement constraints. The 
identification of this type of singularity by mere inspection may be difficult. 


Table III. Five examples that failed to reach optimum weight versus best feasible design 


Number of active constraints at optimum 


Problem 

number 

Optimization 

method 

Percent over- 
design 

versus best feasible design 
Frequency Stress Displacement 

P2 

OC 

5*6 

1 

5 

2 


vs. SQP 

0*0 

1 

8 

2 

P6 

RG 

83-2 

1 

18 

1 


vs. SQP 

0-0 

1 

21 

1 

P7 

SLP 

12 7 

1 

148 

1 


vs. SQP 

0-0 

1 

168 

1 

P8a 

NPSOL 

3-2 

1 

75 

1 


vs. IMSL 

00 

1 

117 

1 

P8c 

IMSL 

7-7 

1 

88 

1 


vs. SUMT 

0*0 

1 

109 

1 
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Table IV. Problems with active constraints exceeding the number of design variables 


r of design 
iables 

Number of active constraints 
at optimum 

10 

11 

6 

20 

12 

170 

57 

119 

10 

11 

5 

9 

16 

19 

16 

35 

25 

38 

25 

30 

49 

75 

49 

76 

24 

28 

24 

28 


Problem number, description 


P2, Tapered 10-bar truss 

P5, Antenna tower 

P7, Geodesic dome 

P8a, Intermediate complexity wing 

P9, 10-bar truss 

P10, 10-bar truss 

PI 7, Cantilever membrane 

P18, Cantilever membrane 

PI 9, 60-bar trussed ring 

P22, 60-bar trussed ring 

P24, Stiffened ring 

P27, Stiffened 60-bar trussed ring 

P30, Stiffened ring 

P33, Stiffened ring 


It is suggested that code developers should address the singularity issue. Singularity alleviation as 
discussed in References 8, 41 and 42 can reduce computation and improve reliability of optimizers. 


Default optimization parameters 

Default parameters (such as convergence criteria, step length, stopping criteria, active con- 
straint region, iteration limitations, etc.) specified by individual optimization codes were used to 
solve the problems. When a problem failed, the default parameters were changed according to the 
instructions specified in the user’s manual of individual codes in an attempt to successfully solve 
the problem. In the solution of the 41 test bed problems, it was necessary to change the default 
optimization parameters quite often in order to reach the correct solution. On an overal basis, 
default parameters of SUMT, SLP, FD, SQP and IMSL algorithms were adequate for the 
solution of most problems. Most of the default parameters for RG and NPSOL were changed to 
improve their performances. 


CONCLUSIONS 

None of the eight optimizers could successfully solve all the problems. Most optimizers, however, 
can solve at least one third of the examples. For large problems, the Cray-YMP CPU time was 
comparable among the optimizers that succeeded. Alleviation of singularity that can occur in 
structural optimization can improve the optimizer efficiency. ... ^ „ 

A single winner which can be called most reliable and efficient could not be identified. Overall, 
three optimizers (IMSL, SUMT and SQP of IDESIGN) scored high marks. For small problems, 
five optimizers (IMSL, SUMT, SQP of IDESIGN, NPSOL and SLP) satisfactorily solved more 
than fifty per cent of the problems. For medium problems, six optimizers (IMSL, SQP o 
IDFSIGN SUMT NPSOL, SLP and FD) produced correct solutions for at least half ol the 
p?offiems For large problems, three optimizers (IMSL, SUMT and SQP of IDESIGN) were 
found to be reliable and efficient. 
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